Number-conserving rate equation for sympathetic cooling of a boson gas 
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We derive a particle number-conserving rate equation for the ground state and for the elementary 
excitations of a bosonic system which is in contact with a gas of a different species (sympathetic 
cooling). We use the Giradeau-Arnowitt method and the model derived by Lewenstein et. al. with 
an additional assumption: the high-excited levels thermalize much faster with the cooling agent 
than the other levels. Evaporation of particles, known to be important in the initial stages of the 
cooling process, is explicitly included. 
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I. INTRODUCTION 

The development of cooling techniques has opened up 
the possibility of studying ultracold gases. In particular, 
the quantum degeneracy of bosons and fermions has been 
investigated. Basically there are two different processes 
which can be used to cool an atomic or molecular ensem- 
ble, evaporative cooling and sympathetic cooling. How- 
ever, there are gases in which the interaction is too weak 
for evaporative cooling to work. In addition, this method 
fails for identical fermions at low temperatures because 
of the exclusion principle Q. In this case, sympathetic 
cooling can be used as an alternative to evaporative cool- 
ing: a cold gas is in thermal contact with another gas to 
be cooled. The first application of the sympathetic cool- 
ing method was in the context of the cooling of charged 
particles Q and later it was extended to neutral atoms 
j|. Only recently has it been applied in the context of 
ultracold atoms. In particular, quantum degeneracy of 
bosons and fermions has been achieved by thermalization 
of atoms of the same species in different internal states 
j| , between two isotopes of the same species [|J , and fi- 
nally between atoms of different species Moreover, 
in contrast to evaporative cooling, sympathetic cooling 
does not lead to a significant loss of either the cooled gas 
or the cooling agent. 

From the theoretical point of view, classical and quan- 
tum models have been derived to describe the dynamics 
of thermalization of one-species and two-species systems. 
In the classical Boltzmann regime, an analytical formula 
for the evolution of the temperature in a mixture of non- 
equal mass atoms and equal mass atoms has been de- 
rived in Refs. [j? 10 1. A quantum mechanical description 
is given in terms of Boltzmann equations |Tl| ] and master 
equations Jl^l5[ . In Ref. (Il]] , sympathetic cooling of an 
atomic Fermi gas by a Bose gas is studied. In particular, 
the equilibrium temperature and the relaxation dynam- 



ics are obtained. In Ref. for the first time, a master 
equation describing the sympathetic cooling of a system 
A is derived, treating the cooling agent B as a heat bath 
at a fixed temperature. More recently, by means of de- 
coherence and ergodicity arguments a rate equation for 
sympathetic cooling of harmonically trapped bosons or 
fermions is obtained ]l3| . In Refs. 0,^5| a detailed quan- 
tum kinetic master equation, which couples the dynamics 
of a trapped condensate to the vapor of noncondensate 
particles, is obtained along with information about the 
growth of a Bose-Einstein condensate. In Refs. [[Hi] pi[ , 
one can find more information about the growth of the 
condensate. 

The purpose of this work is to describe the sympa- 
thetic cooling process of a total-number-conserving sys- 
tem in terms of particles in the ground state and elemen- 
tary excitations. We use the particlc-number-conserving 
Girardeau- Arnowitt formalism ]|a] and the sympathetic 
cooling model developed by M. Lewenstein et al. [ 12) with 
the additional assumption, that, the highly excited lev- 
els of the trap thermalize much faster with the cooling 
agent than the other levels (or in other words, the highly 
excited levels of the trap are in thermal equilibrium with 
the cooling agent). The advantage of the Girardeau- 
Arnowitt formalism is that it covers all cases ranging 
from a total absence of population in the ground state, 
no = 0, to a highly populated ground state no — N. The 
assumption of the fast thermalization of the highly ex- 
cited levels of the trap is related to the separation of the 
levels of the trap into two different bands, in the same 
way proposed by C.W. Gardiner et al. in Refs. n,tt5[, 
the condensate (Be) band and the noncondensate band 
(Bnc), where the latter is in thermal equilibrium with 
the cooling agent. The condensate band includes all trap 
levels which are directly influenced by the presence of the 
condensate. The noncondensate band contains energy 
levels that are sufficiently high for the interaction with 
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the condensate to be negligible. As we will show, the 
master equation for sympathetic cooling can be written 
in terms of two processes: (a) Creation and annihilation 
of particles in the ground state and (b) Creation and 
annihilation of quasiparticles. From this master equa- 
tion we derive a rate equation for this number-conserving 
system under sympathetic cooling which includes the 
growth, scattering and loss (sympathetic evaporation) 
processes. This result is very important, specially now 
that sympathetic cooling is used as an alternative cooling 
method. In addition, this rate equation is more general 
than the one obtained by Gardiner et. al. in refs. p^ , p^| , 
which is only valid in the regime of a highly populated 
condensate [^6| . 

The paper is organized as follows: In section || the de- 
scriptio n o f the sympathetic cooling model is given. In 
section III the master equation describing the dynamics 



of a number-conserving sy stem i n the cond ensat e band 
is derived. In subsections IVA, IV B, and IVC a rate 



equation describing the population growth, the effects of 
the scattered particles as well as the effects of the sym- 
pathetic evaporation (trap losses) is derived. Section |v] 
contains the conclusion. 



II. DESCRIPTION OF THE SYSTEM 

For the sake of completeness and notation we briefly 
review the main results of Refs. jl^Jl^]. However, we 
introduce some modifications to the model presented in 
Rcf . Q : the cooling agent is a gas trapped in a harmonic 
potential, instead of a free gas, and the gas to be cooled 
is confined in an open trap. This situation is closer to 
the actual experimental conditions. 

A system A of Na bosons is subject to sympathetic 
cooling due to its interaction with system B. The cooling 
agent of Ng atoms is in thermal equilibrium at tempera- 
ture Tb- The single-particle states of system B are then 
described by harmonic oscillator wave functions with en- 
ergy eg = Tiujb(£x + (y where u>b is the trap fre- 
quency and £ = {(-xAyAz) are the quantum numbers 
related to the bath. The corresponding Hamiltonian can 
be written 



Hr 



t 



The associated creation and annihilation operators are 
denoted by at and aft, respectively. If, on the other hand, 
the energy is much larger than the trap depth, the system 
A can be approximately described by free particles . In 
this case the states are described by plane waves with en- 
ergy e(R) = h 2 K 2 /(2m), with k denoting the wavevector 
of the particle. 



Ha = ^2 u ™ a n a ™ 



(1) 



where 



U a = 



efi if e <C e t 
c(k) if e > e t 



when e w et, Un cannot be determined in a simple way. 

The two-body interaction between the bath and the 
particles can be written in the form 



V A -b = 



7^ -. 7 .-jata -,bj)i 



For simplicity, we assume a point-like interaction char- 
acterized by a scattering length a. We are interested in 
scattering processes between the bath and tightly bound 
particles, and processes in which low-energy particles of 
the system are scattered into the continuum through the 
interaction with the bath. In the first case, the matrix 
elements of the two-body interaction have the form 



4irh 2 a 
%PAB 



dx^PU^nAxW^-ippix). (2) 



Here, pab = fnxmB I '(m-A + tub) denotes the reduced 
mass. For scattering processes into the continuum, the 
matrix element can be written as 



4irh 2 a 



2mab(2tt) 



3/2 



(3) 



As derived in Ref. [|l2j| . the master equation for the 
time dependence of the reduced density matrix PA(t), ob- 
tained by tracing the total density matrix over the bath, 
has the form 



dp A (t) 
dt 



H A + V A - A ,PA{t) 



+ C PA (t), (4) 



with the creation and annihilation operators for the bath 
&t and bp respectively. 

System A is assumed to be confined in an open trap. 
For energies e much smaller than the trap depth e t , the 
trap can be approximated by a harmonic potential with 
trap frequency uja- The single-particle states in the har- 
monic oscillator potential have the quantum numbers 
n = (n x ,n y ,n z ). The single-particle eigenfunctions are 
labelled ipa(x), the eigenvalues en = Tiu>A{n x + n y + n z ). 



where V A _ A describes the mutual interaction of the 
atoms in system A. The interaction includes the two- 
body interaction Va-a and, in addition, terms corre- 
sponding to shifts produced by the elimination of the 
bath in the master equation. It is assumed that the 
interaction V A - A is weak compared to Va-s, i.e. the 
thermalization of the system is solely determined by the 
interaction with the bath. The equation above is ob- 
tained under the the following assumptions: (i) Markov: 
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the correlation time for the interaction between systems 
A and B is much shorter than the cooling time, and (ii) 
rotating-wave approximation: terms rotating at multi- 
ples of the trap frequency are neglected. 

It was shown in Ref. |l3) that if the number Nb of 
particles in the bath is very large, decoherence acts very 
quickly compared to the equilibration time and reduces 
the density matrix to diagonal form. If only scattering 
processes involving tightly bound atoms of system A are 
involved, the Liouvillian in Eq. (0) becomes 

£PA = X! rm '" 2a !n a nWW4 ,1 'n - a^a^at a S/ cu (i) 

rh^n 



with the rate coefficients 



PA(t)ala^alaft , (5) 



—2 r drV 
xnJrv + 1] exp [i (Aca + Ae B )r//i] (6) 



where n^* is the average occupation number of the heat 
bath oscillator. The energy difference for the system and 
the bath is given by Ae^ — en — em and Ae^ = £/ — en > 
respectively. If particles are scattered into the contin- 
uum, the rate coefficients have similar form but with the 
matrix element defined in Eq. (|J). More details about 
the evaluation of the rate coefficients are given in Refs. 



The condensate band Be includes the ground state 
and all excited levels in the trap which are directly in- 
fluenced by the presence of the condensate. The non- 
condensate band -Bjvc is composed by the highly excited 
levels, which include also the continuum states for e > ej. 
The following assumptions are used in order to derive a 
master equation for the condensate band from Eq. (|^): 

(i) We neglect the correlations between the particles in 

the condensate band and the particles in the non- 
condensate band. Therefore, we assume that the 
complete density matrix of the system A can be 
written as a direct product between the density of 
the condensate band and the density of the non- 
condensate band, i.e., pa — Pc ® Pnc- Since the 
main interest is in the dynamics of the condensate 
band, the noncondensate band is eliminated from 
the description by tracing out the noncondensate 
variable, i.e., pc = Tr NC (pA)- 

(ii) The density matrix of the noncondensate band p^c 

is considered in thermal equilibrium with the cool- 
ing agent at temperature Tb- The levels of this 
band thermalize much faster with the cooling agent 
than the levels of the condensate band. Due to this 
assumption, we will completely neglect collisions 
between particles of the noncondensate band and 
particles of the cooling agent. 

Applying these assumptions to Eq. (||) , a master equation 
for the condensate band can be derived, 



III. MASTER EQUATION FOR THE 
CONDENSATE BAND 



In this section, a master equation describing the dy- 
namics of the condensate band is derived. Basically, 
we extend the formalism developed by Gardiner et al. 
to the case of two distinct gases. The trap-levels 
of the system A are grouped in two bands: the conden- 
sate band (Be) and the band of noncondensed particles 
(Bnc), where the latter is in thermal equilibrium with 
the cooling agent system B (see Figure |l|) . 



trap depth 











Noncondensate Band 


V / 



Condensate Band 



FIG. 1. Distribution of the trap levels in the condensate 
Be and noncondensate B N c bands. 



dpcjt) 
dt 



= Yl r "'™ (N m ) no (2a ] R p c {t)an 
m G B NC 
ft £ B c 



-a n alp c (t) - pc(t)aftal 



(J) 



+ Yl r ™' H <^ + 1 >n (2^p C (*)4 
fn S B NC 
n e B c 



-aianp c (t) - pc(t)aian 



Y I**-* ( 2a^a n pc{t)ala m 
m n 
n, fh £ Be 



(8) 



(9) 
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m ^ ft 
ft,fn £ Be 



-a^an<4ajnPc(*) - PcC*)^ " ^ 01 "^ (10) 



where r™'™ = r^V 3 ^" 6 ") and (N A ) no = 

Tr(ataaPNc)n is the expectation value of the parti- 
cle number operator for state m, given that there are 
no bosons in the ground state. Since the terms with 
fn = ft cancel, this case is already excluded from the 
summation. The first two terms, Eqs.^J^), account for 
collisions between particles in the condensate band and 
in the noncondensate band, and are illustrated in Fig||. 
The loss of particles from the trap, due to their inter- 
action with the cooling agent atoms are also included in 
Eq.(||). The existence of such a process is related to the 
fact that the trap is open, hence the atoms with higher 
energy than the trap energy can indeed escape. In the 
present model these particles are considered to escape 
to continuum states (non-bound states) . In this case the 
single-particle states are described by plane waves. Evap- 
oration is illustrated in Fig.||. The following two terms, 
Eqs.(| 10), describe collisions only between particles in 
the condensate band, which are shown in Fig.^J. 



trap depth 




Noncondensate Band 



Condensate Band 



FIG. 2. Collisions between particles in the condensate band 
Be and in the noncondensate band Bmc 



trap depth 




Noncondensate Band 



Condensate Band 



FIG. 3. Loss of particles from the trap, due to their inter- 
action with the cooling agent atoms. 



trap depth 




Noncondensate Band 



Condensate Band 



FIG. 
B c . 



4. Collisions between particles in the condensate band 



Eqs. (|7|-p~0[) give an accurate treatment of the inter- 
nal dynamics of the condensate band, describing how the 
condensate forms and how the particles are scattered by 
collisions. Effects of elementary excitations of the con- 
densate, representing the "thermal cloud" , do not appear 
explicitly, since they are hidden in the single-particle op- 
erators and at for ft £ Be- In the limit of a highly 
condensed ground state, i.e. no — N, the effects of the 
condensate excitations can be totally neglected |H| . Dur- 
ing the process of condensate formation, however, their 
effects play a crucial role. 

To analyse the role of elementary excitations, we re- 
sort to a formalism of Girardeau and Arnowitt p5[ , 
which conserves the number of particles. The method 
is based on the annihilation fio and creation /3q oper- 
ators of one particle in the ground state, defined as 

A) = (No + l) _1/2 ao and $ = aJ(JV + if 1/2 , and in 
the phonon operators fia, fit defined as fin — fi^aa and 
fit = at fio, where No is the number operator of particles 
in the ground-state, a^ and at are usual annihilation and 
creation operators for the trap level ft. These operators 
obey the Bose commutation relations, 0^, fit,] = 6~ 
The phonons operators can be written in terms of the 
creation and annihilation quasiparticles operators: bft 
and St , i.e., fin = uabn + v a bt and fit = u* n bt + vtb n , 
with \uft\ 2 — |w,t| 2 = 1. The formalism of Girardeau and 
Arnowitt is more general than the one developed by Gar- 
diner |pq| , since it covers all cases ranging from a total 
absence of population in the ground state no = to a 
highly populated ground state no — N. A detailed com- 
parison of both formalisms can be found in Ref . . 

We replace the operators ao, aj, an, at in Eqs. ([7-10) by 
the particle-number-conserving Girardeau- Arnowitt op- 
erators fio , fio ; fin i $ ■ Then, the action of the opertators 
on the density matrix pc(t) is computed. Since the deco- 
herence time is very fast (see discussion in Ref. [|l3| ) , only 
the diagonal terms of the reduced density matrix pc con- 
tribute to the dynamics. Hence, only diagonal terms will 
be considered. In this way, we obtain a master equation 
for the diagonal elements P(r\Q,n) = (no, n\pc(t)\r\Q, ft), 
with ft = (ni, n m , n„, ...) denoting the set of occupa- 
tion number for all quasiparticles: 
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^P(n ,n,f) = ^P g rowth(n ,n,i) + ^P sca tt(n ,n, 



d 
dt 



IjPevapN,™, t). 



(11) 



The rate dP gr owth(no, ft, t)/dt represents the contribution 
of the condensate itself and the quasiparticles for the 
growth of the condensate. The number of particles in the 
condensate band changes. dP sca ±t{r\Q,n,t) / dt describes 
scattered particles in the condensate band. During the 
scattering process the number of particles in the conden- 
sate band does not change. Finally, dP evap ((r\o,n,t) /dt 
includes the evaporation of particles with energies larger 
than the trap depth. 



IV. RATE EQUATIONS FOR THE CONDENSATE 
BAND 

A. Growth 

Consider the diagonal elements of the Eqs. (0— 1|) for 
n = and n ^ with the replacements described in the 
previous section, we arrive at a master equation for the 
occupation probability of the condensate band 



-2^n s r+-(n )P(n ,n) 

ri 

+2 ^2 n « r ^ + ( n o + 1, ft - e! R )P(n Q + 1, ft - e n ) 

ft 

-2 ^2(n H + l)r s + (n , n)P(n , ft), (14) 



with e* n = (...,0,0,1,0,0,...). The coefficients 
r H + ( n o,"0, r «~( n o,n), rt~(n ,7?) and r fl + (n ,?i) are 
defined as 



ri + (n , 


f£) = 


£ r"^(^)n kn| 2 , 








r s ~( n o; 


tf) = 










r±-(n 0: 


ft) = 


E r "'" l W™)n l^| 2 , 






rh^n 


r^ + (no, 


ft) = 


Er™<">,r, + l) n >n| 2 : 



dt 



-P g rowth(no, ft, t) = Rc + Rq, 



(12) 



where Rc represents the growth rate of the condensate 
itself 

R c = 2n r+(n - l)P(n - l,n) - 2(n + l)r+(n )P(n , ft) 
+2(n + l)r-(n + l)P(n + l,ft)- 2n r-(n )P(n , ft). 

(13) 

The coefficients r + (no) and r~(no) carry information 
about particles entering and leaving the condensate 

r+(n )^E r ° Vl ^>no> 
r-(no)^E r ^°>™ + 1 )no- 

m^O 

The rate Rq takes into account the contribution of the 
elementary excitations 

R Q = 2^n s rt + (n - I, ft- e s )P(n - I, ft- en) 

n 

-2£(n« + l)r++(no,n)P(n ,n) 

n 

+2 £(n fl + l)r s "(n + 1, ft + en)P(n Q + 1, ft + en) 

n 

-2Ensrz~P(n ,n) 

n 

+2EKi + l)rt-(n - l,n + e«)P(n - l,n + e«) 



where rh G Bjvc and ft £ Pc- These coefficients describe 
processes of creation and annihilation of phonons inside 
the condensate. 

In order to obtain the dynamics of the population in 
the condensate band, the mean number of particles is 
evaluated. We will start with the contribution of the 
condensate mode itself. In Eq. (^) we use the thermody- 
namic relation between the rate coefficients 

r+(n ) = r-(n )e' 3(eo -' io) 

in combination with the factorization assumption 
P(n , ft) = P(n ) O P(n) = P(n ) . . . ® P(n„) ® . . .. Here 
(3 = 1/}cbTb, eo is the energy of the ground state and 
fiQ = /i(no) the chemical potential for no particles in the 
ground state. For the rate of change of the condensate 
number one then finds 



dt 



| g rowth= 2 ^(n + l)r+(n ) - n r (n ) 

= 2r+(n ) (n (l - e-^«-^) + lj . (15) 



In the limit of no ~ N, with N 3> 1, the equation above 
can be approximated by 



dN 
~~dt 



2T+{N)N[{l~e^ a -^)^, (16) 



where the coefficient T + (N) contains information about 
the sympathetic cooling process between system A and 
cooling agent B. 

If we apply the formalism described here to a gas of 
only one species, we find a rate equation for the growth 
of the condensate which has a form similar to the one 
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obtained by Gardiner and coworkers based on [g6j (see 
Eq. (23) in Ref. g§ and Eq.(168) in Ref. fb§). However, 
due to the properties of the Giradeau-Arnowitt formal- 
ism, the rate equation in our case is not limited only to 
highly populated condensates but also includes conden- 
sates of arbitrary particle number. 

We now explore the contribution of the elementary ex- 
citations for the growth process. We use the thermody- 
namic relation between the coefficients 

r++(n ,n) = r s -(n ,n) e ^-^), 

r+-(no,n) = r-+(no ) n)e^-^), 

and factorize all correlations i.e., P(no,n) = P(no) <8> 
...P(n n )... in Eq. (|l4|) . Thus, the rate equation for the 
elementary excitations is given by 

^Uwth= 2^r++(n 0) n)((l - e^-^K + l) 

-2 £ r+-(n , n) ((I - e^°-^)n n + lj . (17) 

n 

A simplified form of Eq. (Il7|) can be derived by using the 
relation r±"(n , n) = Tt^(n , ft) - rt(n , ft), 

^lgrowth= 2^r+(n ,ri) ((1 - e^-^n* + l) , (18) 

with r+(n 0! n) = E^r«.*<Ar.) no . 

The total growth equation can be obtained from the 
contribution of the condensate mode itself Eq.(|l5|) and 
the elementary excitations Eq.(|l8|). 

B. Scattering 

In this subsection the effects of the scattering processes 
of atoms inside the condensate band are considered, as 
described by Eqs. (|9|— |l0|). During the scattering process 
the number of particles in the condensate band does not 
change. 

We use the same analytical procedure described in 
the previous section to derive a rate equation related to 
the scattering process. However, we will omit details 
of this calculation, since the expressions can be rather 
large. Thus, the corresponding scattering rate equation 
for x\f L — (atcift) has the form 

^Ut= 2 £ n+n [e-^-°)) - n n (l - e^-a))] 

ftjtO 

+2^r^+noMl-e-^- e °)) + l] 



+2^E r^-Mna + i) 

-n s (n, r , + l)e-^~^] 

+2 E r ^T[- n «( n ™ + 1 ) 

+ na (n H + l) e H»<«*-*0] 

+2 £ r-++[(n^ + l)(n^ + l) 

_ n . n . p -fc^)l 

+ (n„ 1 + l)(n K + l)e-' 3 ^-^], (19) 

with 

r«+(n ,n) = r°^| 2 , 
r^(n 0) n) = r^K| 2 , 

t^uu-\ — /„ ->\ pn.mi |2i i2 

r n,^ ( n 0:^) = T M |«rft| , 

tiwi'H — /„ c*\ -pin, mi |2i i2 

r n,^ ("0,"-) = r |W S \Vfh\ , 

r S + ("o,n) =r"^|^| 2 |^[ 2 , 
r!;+ + (n 0l «)Er^|, fi | 2 |^| 2 (20) 

The first two terms of Eq. ( |l9| ) account for the scattering 
between the particles in the ground state and the quasi- 
particles of the n-level. The other terms describe the 
scattering between quasiparticles of different levels. 

C. Evaporation 

We now discuss the evaporation of particles from the 
trap. Atoms can escape from the condensate band of 
the trap to unbound states. This loss of particles from 
the trap is induced by their interaction with the cool- 
ing agent. This process is particularly important in the 
initial stage of the cooling process Jt],[§. To implement 
evaporation in the present model, particles are assumed 
to escape to continuum states (unbound states) inside 
the noncondensate band. Only Eq. (g) contributes to 
this process, since the rate coefficients of Eq. (^), which 
includes the overlap of the wave functions involved in the 
process, is completely negligible. Taking the diagonal el- 
ements of Eq. (||), Povap is obtained 

^Pova P (no,«,i) = -2n 7~(no)-P(n ,™) 
+2(n + l) 7 -(n + l)P(n + l,n) 

+2 E {nn + lh n ~(no + l,n + e n )P{n + l,n + e n ) 
n e B c 

-2 n ? -j7-~(n )P(n ,n) 
n e B c 



G 



-2 ^2 n «7 H + ( n o + l,n- e n )P(n + l,n- e n ) 

u e b c 

-2 J] (n t1 + l) 7 -+(no,n)P(n ,n), (21) 
ra G B c 



where 



-(n ,n)= ^ P 



meB 



(t) 

NC 



7 r-(n ,n)^ £ r"^| UH | 2 , 



(t) 

NC 



7 r+(n ,n)= J] r ™'"l^| 2 

™6 B NC 

= 7« + (no,n) -7 S 



(22) 



with B^ c denoting the noncondensate band with energy 
e-m > e t- As mentioned at the end of section 55, these rate 
coefBcients r m >° and T m,n have a different form from the 
ones which appear in the growth and scattering processes, 
since now m describes unbound states. The first line of 
Eq.(pH) describes the evaporation of an atom out of the 
ground state while the others describe the evaporation 
of an atom out of the n-excited level of the condensate 
band. 

Following the same procedure described in the previ- 
ous subsections, a rate equation for the evaporation of 
particles from the ground state and the evaporation of 
quasiparticles in the condensate band is obtained: 



dt |cvap : 

dt lovap 



-2 7 (n ,n) n | cvap , 
- 2 H7n(no,n) n«| evap , 



with 



7a 



X 1 pTTl.n 

/ j n,fn' 
m£B« 



(23) 



(24) 



for the thermalization of a system A which is in thermal 
contact with a bath B. All the information about the 
dynamics of the thermalization process is contained in 
the coefficients T. Analytical and numerical evaluation 
of the coefficients for a specific system can be carried out 
along the lines described in Refs. Jl2] , ^3[ . 

The result obtained here can be applied to describe 
the dynamics of sympathetic cooling of a gas in thermal 
contact with a cooling agent, in terms of the population 
in the ground-state and elementary excitations. In par- 
ticular, it can be aplied to the case where the cooling 
agent thermalize much faster compared to the thermal- 
ization of the system. The description remains valid for 
the quantum degenerate regime T <T C , where T c is the 
critical temperature of the gas. 

The use of the Giradeau-Arnowitt method has opened 
up the possibility to describe sympathetic cooling in 
terms of particles in the ground state and elementary 
excitation during the whole cooling process, i.e., from 
the case of total absence of particles in the condensate 
(no = 0) to the case of a highly populated ground state 
(no = N). In addition, we include the effects of evapo- 
ration which are especially important during the initial 
stages of the cooling process. To our knowledge this is 
the first time that a formal complete number-conserving 
description of sympathetic cooling in terms of particles 
in the ground state and elementary excitations is given. 
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V. CONCLUSION 

The total rate equation describing a number- 
conserving population of the condensate band can be ob- 
tained from Eqs.rfjM^HIl^ 



dnn 
~~dl 



dn R 



dn R 



dux , 



^ + ^ (growth ~i~ ^ + Iscatt ~t~ ^ + levapi (^5) 



dt 



dt 



where no and n^ denote, respectively, the population in 
the condensate and the population of the elementary ex- 
citations, which are commonly referred to as "thermal 
cloud" . The equation above gives a complete description 
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